Microscopic mechanisms of spin-dependent electric polarization in 3d oxides 
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We address a systematic microscopic theory of spin-dependent electric polarization in 3d oxides 
starting with a generic three-site two-hole cluster. A perturbation scheme realistic for 3d oxides 
is applied which implies the quenching of orbital moments by low-symmetry crystal field, strong 
intra-atomic correlations, the dp-transfer effects, and rather small spin-orbital coupling. An ef- 
fective spin operator of the electric dipole moment is deduced incorporating both nonrelativistic 
tx (§i • §2) and relativistic oc [si x sa] terms. The nonrelativistic electronic polarization mechanism 
related with the effects of the redistribution of the local on-site charge density due to pd covalency 
and exchange coupling is believed to govern the multiferroic behaviour in 3d oxides. The relativis- 
tic exchange- dipole moment is mainly stems from the nonrelativistic one due to the perturbation 
effect of Dzyaloshinsky-Moriya coupling and is estimated to be a weak contributor to the electric 
polarization observed in the most of 3d multiferroics. 



I. INTRODUCTION 

Strong coupling of magnetism and ferroelectricity was 
recently uncovered in rare earth manganites with the gen- 
eral formula RMnC>3 and RM^Os, where R= a rare 
earth ion, or Y (see e.g., RefsU and review articles 
RefsHfl!)- In magnetically ordered state below these 
ferroelectric magnets, or multiferroics, exhibit an excep- 
tionally strong sensitivity to an applied magnetic field, 
which induces reversals and sudden flops of the electric 
polarization vector, and results in a strong enhancement 
of dielectric constant. Vice versa also an applied electric 
field affects the magnetic properties such as the helicity. 

Since the Astrov's discovery of magnetoelectric effect 
in C^O^ there were proposed several microscopic mech- 
anisms of magnetoelectric coupling^ however, the mul- 
tiferroicity have generated an impressive revival of the 
activity in this field. Currently two essentially different 
spin structures of net electric polarization in crystals are 
considered: i) a bilinear nonrelativistic symmetric spin 
coupling££i£ 

p s = ^n^„(s m -s„) (i) 

mn 

or ii) a bilinear relativistic antisymmetric spin 
coupling 9 - ' 10 ! 11 

Pa = J2 I 11 ™ X t S " X S «l] > ( 2 ) 

mn 

respectively. The effective dipole moments 11^* depend 
on the m, n orbital states and the mn bonding geometry. 

If the first term stems somehow or other from a 
spin isotropic Heisenberg exchange interaction (see, e.g. 
Refslall 2|). the second term does from antisymmetric 
Dzyaloshinsky-Moriya (DM) coupling. Namely the sec- 
ond, or "spin-current" term is at present frequently 
considered to be one of the leading mechanisms of 
multiferroicity i 9 ' 10 ! 11 ! 13 ! 14 ! 15 ! 16 However, there are no- 
table exceptions, in particular the manganites RM112O5, 



H0M11O3, where a ferroelectric polarization can appear 
without any indication of magnetic chiral symmetry 
breaking^, and delafossite CuFei_^Al :r 02, where the 
helimagnetic ordering generates a spontaneous electric 
polarization || to the helical axisil, in sharp contrast with 
the prediction of the spin current model. 

The recent observations of multiferroic behaviour con- 
comitant the incommensurate spin spiral ordering in 
chain cuprates LiCuVC>4 by Naito et al^ and LiCu202 
by Park et al^ challenge the multiferroic community. 
At first sight, these quantum s=l/2 ID helicoidal mag- 
nets seem to be prototypical examples of ID spiral- 
magnetic ferroelectrics revealing the relativistic mech- 
anism of "ferroelectricity caused by spin-currents". 9 In- 
deed, Naito et al*ZL claim that the electric polarization 
in LiCuV04 can be understood by the relation predicted 
by spin-current models^^ However, LiCu202 shows up 
a behavior which is obviously counterintuitive within 
the framework of spiral-magnetic ferroelectricitypiS More- 
over, in contrast with Park et al^, Naito et al^ have 
not found any evidence for a ferroelectric transition in 
LiCu2C>2. 

The ferroelectric anomaly in LiVCu04 reveals a mag- 
nitude comparable or even larger than that of the multi- 
ferroic M3V2O8 where the magnetic ordering drives the 
electric polarization w 10 2 fiC/m 2 (Refl2fJ) that rep- 
resents a typical magnitude of polarization induced by 
magnetic reordering in multiferroics. However, such an 
anomalously strong magnetoelectric effect seems to be 
an unexpexted feature for a system with e g -holes and 
a nearly perfect highly symmetric chain structure with 
the edge-shared CUO4 plaquettes which both are un- 
favourable for a strong spin-electric coupling. Thus the 
giant magnetoelectric effect in the title cuprate raises fun- 
damental questions about its microscopic origin. 

Microscopic quantum theory of ME effect has not yet 
been fully developed, although several scenarios for par- 
ticular materials have been proposed based on the effec- 
tive spin HamiltonianiiiSiiii In a recent paper Katsura et 
al£ presented a mechanism of the giant ME effect the- 
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oretically derived "in terms of a microscopic electronic 
model for noncollinear magnets" . The authors derived 
the expression for the electric dipole moment for the spin 
pair as follows: 

Pij = a [Rij x [Si x Sj}} , (3) 

where Rjj denotes the vector connecting the two sites i 
and j, Sij are spin moments, a is an exchange-relativistic 
parameter. It is worth noting that the mechanism 
also implies the uniform polarization accompanying the 
spin-density wave. However, the original "spin-current" 
model by Katsura et al£ and its later version s 14 ' 15 ' 16 
seem to be questionable as the authors proceed with an 
unrealistic scenario. Indeed, when addressing a generic 
M1-O-M2 system they groundlessly assume an effective 
Zeeman field to align noncollinearly the spins of 3d- 
electrons and to provide a nonzero value of the two-site 
spin current [Si x S2]. To justify their approach, the 
author s 9 ' 14 ! 15 ' 16 are forced to assume a colossal magni- 
tude of this fictious field resulting in an enormous Zeeman 
splitting of several eV. Second, Katsura et aL— start with 
an unrealistic for 3d-oxides strong spin-orbital coupling 
limit for ti g electrons 21 - which formally implies A 3> U 
and a full neglect of the low-symmetry crystal field and 
orbital quenching effect.— The authors 9 - ' 14 ' 15 ' 16 do heav- 
ily (up to two orders of magnitude!) overestimate the 
numerical value of the overlap dipole matrix element 
/(Rdp) = f d yz (r)yp z (r + H,dp)dr which defines maximal 
value of respective electric dipole moments. It seems, the 
authors ignore the well developed techniques to proceed 
with pd-covalency, exchange, and spin-orbital coupling in 
3d oxides. 

Alternative mechanism of giant magnetoelectricity 
based on the antisymmetric DM type magnetoelastic cou- 
pling was proposed recently by Sergienko and Dagotto. 11 
However, here we meet with a "weak" contributor. In- 
deed, the minimal value of 7 parameter (7 = dD/dll) 
needed to explain experimental phase transition in mul- 
tiferroic manganites is two orders of magnitude larger 
than the reasonable microscopic estimations. 11 

In our opinion, a misunderstanding exists regarding 
the relative role of the off-center ionic displacements (lat- 
tice effects) and electronic contributions to a resultant 
electric polarization. Many authors consider the giant 
multiferroicity requires the existence of sizeable atomic 
displacements and structural distortionsj 24 ' 25 One would 
expect a transition to a structure with polar symmetry to 
occur at the onset of ferroelectricity, but neutron diffrac- 
tion studies thus far have failed to find direct evidence of 
such changes. 26 Earlier synchrotron x-ray studies found 
some evidence of lattice modulation in the ferroelectric 
phase of YMn2 05^21 though the atomic displacements 
seem to be extremely small. Other structural works 
have not reported any signature of atomic displacements 
~ O.OOlA at the ferroelectric phase transition which can 
explain the polarization observed in this family of com- 
pounds. This questions the microscopic model by Harris 
et al£^ supposing the dominant role of the displacement 



derivatives of the exchange integrals, especially because 
the Bloch's rule -fy^ ~ 10 (Refill) point to magni- 
tudes of these derivatives as insufficient to explain the 
~ O.OOlA displacements. However, several phonons in 
TbMn2 05 exhibit clesr correlations to the ferroelectric- 
ity of these materials.— The signatures of the loss of in- 
version symmetry in the ferroelectric phase were found 
by the appearance of a infrared active phonon that was 
only Raman active in the paraelectric phase. A seem- 
ing contradiction, we think a result of an oversimplified 
approach to the lattice dynamics. Indeed, the effects of 
nuclear displacements and electron polarization should 
be described on equal footing, e.g., in frames of the well- 
known shell model of Dick and Overhauser— widely used 
in lattice dynamics. In frames of the model the ionic con- 
figuration with filled electron shells is considered to be 
composed of an outer spherical shell of 2(21+1) electrons 
and a core consisting of the nucleus and the remaining 
electrons. In an electric field the rigid shell retains its 
spherical charge distribution but moves bodily with re- 
spect to the core. The polarizability is made finite by 
a harmonic restoring force of spring constant k which 
acts between the core and shell. The shells of two ions 
repel one another and tend to become displaced with 
respect to the ion cores because of this repulsion. Shell 
and core displacements may be of comparable magnitude. 
The conventional shell model does not take into account 
the spin and orbital degrees of freedom, hence it cannot 
describe the multiferroic effects. In fact, the displace- 
ments of both the atomic core and electron shell would 
depend on the spin surroundings producing the sinergetic 
effect of spin-dependent electric polarization. Obviously, 
this effect manifest itself differently in neutron and x-ray 
diffraction experiments. Sorting out two contributions is 
a key issue in the field. 

The authors of recent papers Refs l3l"ll32l made use 
of first principles calculations to examine the spin- 
dependent electric polarization in the orthorhombic mul- 
tiferroic HoMnOg^ii and in spin spiral chain cuprates 
LiCuV04 and LiCu202.— However, their results are 
highly questionable since the basic starting points of the 
current versions of such spin-polarized approaches as the 
LSDA exclude any possibility to obtain a reliable quanti- 
tative estimation of the spin-dependent electric polariza- 
tion in multiferroics. Indeed, the basic drawback of the 
spin-polarized approaches is that these start with a local 
density functional in the form (see, e.g. Ref[H) 

m(r) 

v(r) = v n(r) + Au n(r), m(r) ]{& ■ - — — -) , 

|m(r)| 

where n(r), m(r) are the electron and spin magnetic den- 
sity, respectively, a is the Pauli matrix, that is these im- 
ply presence of a large fictious local one-electron spin- 
magnetic field cx (vr — v*), where v*'* are the on-site 
LSDA spin-up and spin-down potentials. Magnitude of 
the field is considered to be governed by the intra-atomic 
Hund exchange, while its orientation does by the effec- 
tive molecular, or inter-atomic exchange fields. Despite 
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the supposedly spin nature of the field it produces an 
unphysically giant spin-dependent rearrangement of the 
charge density that cannot be reproduced within any con- 
ventional technique operating with spin Hamiltonians. 
Furthermore, a direct link with the orientation of the 
field makes the effect of the spin configuration onto the 
charge distribution to be unphysically large. Similar ef- 
fects cannot be reproduced in frames of any conventional 
Heisenberg model. In general, the LSDA method to han- 
dle a spin degree of freedom is absolutely incompatible 
with a conventional approach based on the spin Hamil- 
tonian concept. There are some intractable problems 
with a match making between the conventional formal- 
ism of a spin Hamiltonian and LSDA approach to the ex- 
change and exchange-relativistic effects. Visibly plausi- 
ble numerical results for different exchange and exchange- 
relativistic parameters reported in many LSDA investi- 
gations (see, e.g., RefsHjHH) evidence only a potential 
capacity of the LSDA based models for semiquantitative 
estimations, rather than for reliable quantitative data. It 
is worth noting that for all of these "advantageous" in- 
stances the matter concerns the handling of certain clas- 
sical Neel-like spin configurations (ferro-, antiferro-, spi- 
ral,...) and search for a compatibility with a mapping 
made with a conventional quantum spin Hamiltonian. 
It's quite another matter when one addresses the search 
of the charge density redistribution induced by a spin 
configuration. In such a case the straightforward applica- 
tion of the LSDA scheme can lead to an unphysical over- 
estimation of the effects or even to qualitatively incorrect 
results due to an unphysical effect of a breaking of spatial 
symmetry induced by a spin configuration. As an exam- 
ple, we refer to the papers by Picozzi et alM- and Xiang 
and Whangbo^ where the authors made use of the first- 
principle LSDA calculations to study the microscopic ori- 
gin of ferroelectricity induced by magnetic order in or- 
thorhombic HoMnOs and in quasi-lD cuprates LiCu2C>2 
and LiCuV04, respectively. The calculated total nonrel- 
ativistic polarization of the AFM-E phase in HoMnOs 
exceeds the experimentally measured one by more than 
three orders of magnitude. In terms of a conventional 
scheme the AFM-E ordering turns out to be accompanied 
by a colossal exchange striction of the order of several 
percents that exceeds all the thinkable magnitudes (see 
Table I in Refl3ll). The relativistic LSDA calculations for 
the optimized structures of quasi-lD cuprates^ yield the 
results that disagree with experiment both quantitatively 
and qualitatively. Again we see an unphysically strong 
overestimation of the spin-induced electric polarization. 
Interestingly, that the making use of experimental cen- 
trosymmctric structures leads to a strong suppression 
by order of magnitude of the calculated polarizations, 
clearly confirming the unphysically strong LSDA over- 
estimation of spin-induced structural and charge density 
distortions. Summarizing, we should emphasize two weak 
points of so-called first-principle calculations which ap- 
pear as usual to be well forgotten in the literature. First, 
these approaches imply the spin configuration induces 



immediately the appropriate kinematic breaking of spa- 
tial symmetry that makes the symmetry-breaking effect 
of a spin configuration unphysically large. Conventional 
schemes imply just opposite, however, a physically rea- 
sonable picture when the charge and orbital anisotropics 
induce a spin anisotropy. Second, these neglect quan- 
tum fluctuations, that restricts drastically their appli- 
cability to a correct description of the high-order per- 
turbation effects. Overall, the LSDA approach seems 
to be more or less justified for a semiquantitative de- 
scription of exchange coupling effects for materials with 
a classical Neel-like collinear magnetic order. However, 
it can lead to erroneous results for systems and effects 
where the symmetry breaking and quantum fluctuations 
are of a principal importance such as: i) noncollinear 
spin configurations, in particular, quantum s=l/2 mag- 
nets, ii) relativistic effects, such as the symmetric spin 
anisotropy, antisymmetric DM coupling, and, iii) spin- 
dependent electric polarization. Indeed, a correct treat- 
ment of these high-order perturbation effects needs in a 
correct account both of local symmetry and of quantum 
fluctuations(see, e.g., Refl36l). 

It is worth noting that the spin-current scenario by 
Katsura et alM- starts with the same LSDA-like as- 
sumption of unphysically large symmetry-breaking spin- 
magnetic field. Surprisingly, despite the problems with 
the model validation and quantitative estimations the 
spin-current mechanism is currently addressed to be 
responsible for the emergence of ferroelectric polariza- 
tion in new multifcrroics such as orthorhombic RMnC>3, 
Ni 3 V 2 8 , MnW0 4 , CoCr 2 4 , CuFe0 2 where the inver- 
sion symmetry breaking is related to noncollinear spiral 
magnetic structures^ "Ferroelectricity caused by spin- 
currents" has established itself as the leading paradigm 
for both theoretical and experimental investigations in 
the field of strong multiferroic coupling. However, a 
"rule" that chiral symmetry needs to be broken in or- 
der to induce a ferroelectric moment at a magnetic phase 
transition is questionable. Moreover, there are increas- 
ing doubts whether weak exchange-relativistic coupling 
can generate giant electric polarization observed in mul- 
tiferroics. Thus we may assert that a true microscopic 
mechanism of giant magnetoelectric effect is still miss- 
ing. 

Below we propose a systematic microscopic theory of 
spin-dependent electric polarization which implies the 
derivation of effective spin operators for nonrelativistic 
and relativistic contributions to electric polarization of 
the generic three-site two-hole cluster such as Cui-O- 
Cu 2 and does not imply any fictious Zeeman fields to 
align the spins. We make use of conventional well-known 
approaches to account for the p<i-covalent effects, intra- 
atomic correlations, crystal field, and spin-orbital cou- 
pling. Despite the description is focused on Cui-0-Cu2 
clusters typical for different cuprates, the generalization 
of the results on the M1-O-M2 clusters in other 3d oxides 
is trivial. 

The paper is organized as follows. In Sec. II we con- 
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sider the effects of pd covalency and spin-orbital coupling 
in a three-site two-hole Cui-0-Cu 2 cluster. Nonrelativis- 
tic and relativistic mechanisms of spin-dependent electric 
polarization with local and nonlocal terms are discussed 
in Sees. Ill and IV, respectively. In Sec.V we address 
an alternative approach to nonrclativistic mechanism of 
spin-dependent electric polarization induced by a parity- 
breaking exchange interaction. In Sec. VI we show a lack 
of the spin-dependent electric polarization effects for an 
isolated CuC>2 chain. 



II. THREE-SITE TWO-HOLE M1-O-M2 
CLUSTER 

Before proceeding with electric polarization effects we 
address the generic three-site M1-O-M2 cluster which 
forms a basic element of crystalline and electron structure 
for 3d oxides. A realistic perturbation scheme needed to 
describe the active M 3d and O 2p electron states implies 
the strong intra-atomic correlations, the comparable ef- 
fect of crystal field, the quenching of orbital moments by 
low-symmetry crystal field, account for the dp-transfer 
up to the fourth order effects, and rather small spin- 
orbital coupling. To this end we make use of a tech- 
nique suggested in refsj 36 i 37 to derive the expressions 
for the copper and oxygen spin-orbital contributions to 
Dzyaloshinsky-Moriya coupling in copper oxides. For il- 
lustration, below we address a typical for cuprates the 
three-center (CU1-O-CU2) two-hole system with tetrago- 
nal Cu on-site symmetry and ground Cu 3d x 2_ y 2 states 
(see Fig. 1) which conventional bilinear spin Hamiltonian 
is written in terms of copper spins as follows 

H S (12) = J 12 (si • s 2 ) + D12 ■ [sx x s 2 ] + S1K12 s 2 , (4) 

where J12 > is an exchange integral, D12 is the 

Dzyaloshinsky vector, K12 is a symmetric second-rank 
tensor of the anisotropy constants. In contrast with 

</i2,Ki2, the Dzyaloshinsky vector D12 is antisymmet- 
ric with regard to the site permutation: D12 = D21. 




FIG. 1: Geometry of the three-center (Cu-O-Cu) two-hole 
system with ground Cu 3d x 2 _ y i states. 



Hereafter we will denote J12 = J, K12 = K, D12 = D, re- 
spectively. It should be noted that making use of effective 
spin Hamiltonian ^ implies a removal of orbital degree 
of freedom that calls for a caution with DM coupling as, 
strictly speaking, it changes both a spin multiplicity, and 
an orbital state. 

For a composite two s — 1/2 spins system one should 
consider three types of the vector order parameters: 

S = Sl +s 2 ; V = s 1 -s 2 ;T = 2[s 1 xs 2 ] (5) 

with a kinematic constraint: 

S 2 +V 2 = 3i; (S-V) = 0; (f -V) = 6i; [txV] = S. (6) 

Here S is a net spin of the pair, the V operator describes 
the effect of local antiferromagnetic order, or staggered 
spin polarization, while T operator may be associated 
with a pair spin current. Both T and V operators change 
the net spin multiplicity with matrix elements 

(00|f m |ln) = -(ln|f m |00) = i5 mn ; 

(00\V m \ln) = (ln|V> m |00) = S mn , (7) 

where we made use of Cartesian basis for 5 = 1. The 
eigenstates of the operators V n and T„ with nonzero 
eigenvalues ±1 form Neel doublets ^75 (1 00) ± |ln)) and 

DM doublets ^=(|00) ± i|ln)), respectively. The Neel 
doublets correspond to classical collinear antiferromag- 
netic spin configurations, while the DM doublets corre- 
spond to quantum spin configurations which sometimes 
are associated with a rectangular 90° spin ordering in the 
plane orthogonal to the Dzyaloshinsky vector. 

It should be noted that the spin Hamiltonians can be 
reduced to within a constant to a spin operator acting in 
a net spin space 

H s = ij(S 2 -V 2 ) + i(D-T) + isK S S-ivK V V. (8) 

4 2 4 4 

Hereafter we assume a tetragonal symmetry at Cu sites 
with local coordinate systems as shown in Fig[TJ The 
global xyz coordinate system is chosen so as CU1-O-CU2 
plane coincides with xy-plane, x-axis is directed along 
CU1-CU2 bond. In such a case the basic unit vectors 
x, y, z can be written in local systems of Cui and Cu2 
sites as follows: 

x = (sin — , — cos — cosoi, — cos — smoi); 



y = (cos - , sin - cos 5%, sin - sinJi); z = (0, sin<5i, cos Si) 

for Cui, while for CU2 site 0,Si should be replaced by 
—6,62, respectively. 

We start with the construction of spin-singlet and spin- 
triplet wave functions for our three-center two-hole sys- 
tem taking account of the p-d hopping, on-site hole-hole 
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repulsion, and crystal field effects for excited configura- 
tions {n} (Oil, 110, 020, 200, 002) with different hole 
occupation of Cui, O, and Cu2 sites, respectively. The 
p-d hopping for Cu-0 bond implies a conventional Hamil- 
tonian 

Hp d = E tpadfiPadp + h.c. , (9) 

a/3 

where p\ creates a hole in the a state on the oxygen 
site, while dp annihilates a hole in the (3 state on the 
copper site; t pa d(3 is a pd-transfer integral (ip a .d x2 _ B 2 = 

'^~tp z d z 2 — tpda > 0, tpyd^y = tpdn > 0). 

For basic 101 configuration with two d x i_ y i holes lo- 
calized on its parent sites we arrive at a perturbed wave 
function as follows 

*101;SM =VS[®W1;SM+ ^ c {™} ( 25+1 r)$ {n}; rSAf] , 
r{n}^101 

(10) 

where the summation runs both on different configura- 
tions and different orbital T states; 

vs = (i+£lW 2S+1 r)l 2 ) (ii) 

\ {«}r / 

is a normalization factor. It is worth noting that 
the probability amplitudes, or hybridization parameters, 
C{oii},c { iio} °c ipd,C{2oo},c { 020},C{oo2} qc t 2 pd . For in- 
stance, 

c ^=-4 ^fcy sh 4 ; (i2) 
c ^ )= -4e^ c °4' (i3) 

where c Btt {dp) = c uo (dp), c 8)t (pd) = c Q11 (dp) are proba- 
bility amplitudes for different singlet (c s ) and triplet (q) 
110 (Cu\2id x 2_ y 202p Xy y) and 011 {02p Xi yCu2&d x 2_ y 2) 
configurations in the ground state wave function, respec- 
tively; c s ,t(dp x )=-c Stt (p x d), c St t(dp y ) = c„ t t(pyd), t dpa 
is a hole dp-transfer integral. The energies E s ^(dp x ,y) 
are those for singlet and triplet states of dp X:V configura- 
tions, respectively: E s>t (dp Xty ) = e x>y + K dpX} y ± 
where Kd pXlV and Idpx,y are Coulomb and exchange dp- 
integrals, respectively.lt is worth noting that the energies 
e x<y accomodate both the pd transfer energy A p d and 
crystal field effectsie^.y = A p d + 5e Xty . To account for or- 
bital effects for Cui ; 2 3d holes and the covalency induced 
mixing of different orbital states for 101 configuration we 
should introduce an effective exchange Hamiltonian 

He* = \ E J(uPl5)d\ ct ^d\^ l d 2l »di5v,' + h.c. (14) 

a/37<5pp' 

Here d\ all creates a hole in the crth 3d orbital on Cui 
site with spin projection [i. Exchange Hamiltonian (114[) 



involves both spinless and spin-dependent terms, how- 
ever, it preserves the spin multiplicity of CU1-O-CU2 sys- 
tem. Exchange parameters J(af3jS) are of the order of 
tp d . The conventional exchange integral can be written 
as follows: 

J = E [\c { n}( 3 n 2 E 3T ({n})-\c {n} er)\ 2 Ei r ({n})}. 

{n},T 

(15) 

To account for relativistic effects in the three-site clus- 
ter one should incorporate the spin-orbital coupling both 
for 3d- and 2p-holes. Local spin-orbital coupling is taken 
as follows: 

V so = E^' 1 '^ ' S = 



^-[(h +h)-S + (h - i 2 ) ■ V] = A s ■ S + A v ' ■ V (16) 

with a single particle constant £ n ; > for electrons 
and £ n i < for holes. We make use of orbital 
matrix elements: for Cu 3d holes (d x 2_ y 2\l x \d yz ) = 
(d x 2_ y 2\l y \d xz ) = i,(d x 2^ y 2\l z \d xy ) = -2i, (i\lj\k) = 
-itijk with Cu 3d yz = |l), 3d xz =|2) 3d xy =\3), and for O 
2p holes (pi\lj\pk) = itijk- Free cuprous Cu 2+ ion is de- 
scribed by a large spin-orbital coupling with |£ 3t ;| = 0.1 
eV (see, e.g., Refl38l). though its value may be signifi- 
cantly reduced in oxides. Information regarding the £ 2 p 
value for the oxygen O 2- ion in oxides is scant if any. 
Usually one considers the spin-orbital coupling on the 
oxygen to be much smaller than that on the copper, 
and therefore may be neglected. 39 ! 40 However, even for 
a free oxygen atom the electron spin orbital coupling 
turns out to reach of appreciable magnitude: £ 2p — 0.02 
eV^i while for the oxygen 2 ~ ion in oxides one ex- 
pects the visible enhancement of spin-orbital coupling 
due to a larger compactness of 2p wave function^ 2 - If to 
account for oc (r~ 3 ) n i and compare these quantities 
for the copper and the oxygen ((r^ 3 )^ ~ 6 — 8 a.u. and 
( r ~ 3 )2p ~ 4 a.u., respectively^ 2 -) we arrive at a maximum 
factor two difference in £3^ and £?p (see, also Refl43h. 
The Dzyaloshinsky-Moriya coupling 

Hdm = D 12 • [Si x s 2 ] = i(D • f ) (17) 

can be addressed to be a result of a projection of the spin- 
orbital operator Vso = Vso(Cui) + Vso{0) + Vso{Cu 2 ) 
on the ground state singlet-triplet manifold. 36 Remark- 
ably that the net Dzyaloshinsky vector D12 has a par- 
ticularly local structure to be a superposition of partial 
contributions of different ions [i = 1, 0, 2) and ionic con- 
figurations {n} = 101, 110, 011, 200, 020, 002 

D = J2 D l n> • (18) 

i,{n} 

The partial contributions are analyzed in details in 

RefM 



6 



III. NONRELATIVISTIC MECHANISM OF 

SPIN-DEPENDENT ELECTRIC 
POLARIZATION:LOCAL AND NONLOCAL 
TERMS 

Projecting electric dipole moment P = |e|(ri + on 
the spin singlet or triplet ground state of two-hole system 
we arrive at an effective electric polarization of three- 
center system (P) s = (*ioi ; sa/|P|*ioi ; sm) to consist of 
local and nonlocal terms:P = V local + V nonlocal , which 
accomodate the diagonal and nondiagonal on the ionic 
configurations matrix elements, respectively. The local 
contribution describes the redistribution of the local on- 
site charge density and can be written as follows: 

(P) l s ° cal = leH^I 2 [(Ri + R 2 + (Rj + R ) ]T |c 110 (^r)| 2 



-(Ro + R 2 J2 |con(ST)| 2 + 2R G \ c 02o(ST)f 



-2R! |c 2 oo(5T)| 2 + 2R 2 ^ |c 002 (^r)| 2 



(19) 



where Po = |e|(Ri + R 2 ) is a bare purely ionic two-hole 
dipole moment. This dipole moment incorporates both 
the large (cx t pd ) and small (oc t pd ) contributions. Obvi- 
ously, the net local electric polarization can be expressed 
as a sum of local dipole moments: 



(pyocal 



/- D \local 



though, from the other hand, it is easy to show that it 
depends only on R^ vectors (Rio, R207 Ri2)- To this 
end one should carefully proceed with the normalization 
factor in (I19|) . It is worth noting that the net local electric 
polarization lies in the CU1-O-CU2 plane. 

The nonlocal, or overlap contribution is related with 
nondiagonal two-site matrix elements of P and in the 
lowest order with respect to a pd transfer integral can be 
written as follows: 



nonlocal 



2\e\n s 



[c S (p l d)(2p l \r\3d[ 



(i) 



c s (dpi)(2pi\r\3d^ 



(2) 



-i—x,y 



(20) 



(Px)s,t = — HeKcos 2 8 2 - cos 2 <5i)sin( 



(2p y \y\id x 2_ y 2) 



E s j(dp x ) E s j(dp y ) 



(21) 



(P y ) s j = -y/3\e\t 



,„la COS-[( 



cos 2 Si 



'■s 2 ) 



(2p v \y\3d x 



>sin 2 -( 



1 



2 y E Stt (dp x ) E St t{dp y )- 



2(2p x \x\3d x 2 



sm 



'E Sit (dp y ) E s j(dp v 



:)] . ( 22 ) 



where all the matrix elements are taken in local coordi- 
nates of Cu sites. For a symmetric d-orbitals arrange- 
ment with <5i = 82 the a;-component of electric polariza- 
tion {P x ) s ,t turns into zero regardless the bonding angle 
8, whereas the y-component {P y ) St t turns into zero only 
if 6 = 7T, that is for collinear Cu-O-Cu bonding. It should 
be noted that both the partial and net nonlocal contribu- 
tions to electric polarization lie in the CU1-O-CU2 plane 
and are believed to have the same symmetry properties. 

Nominally, the nonlocal contribution to the electric 
dipole moment is proportional to the pd transfer integral, 
however, actually the two-site dipole matrix elements in- 
directly are proportional to the pd overlap integral S p d 
that in a sense equalizes the nonlocal and local terms. 
Let address the problem of the two-site dipole matrix el- 
ements in more details because their correct estimation 
allows to make a reliable conclusion regarding the rela- 
tion between local and nonlocal terms, and the resultant 
effect itself. For instance, Katsura et alM did heavily (up 
to two orders of magnitude!) overestimate the numerical 
value of the integral I(Rd p ) = / d yz (r)yp z (r + Rd p )dr 
which defines maximal value of respective electric dipole 
moments. Indeed, the authors erroneously replaced the 
actually two-site integral by a respective one-site integral 
with the hydrogen-like 3d- and 2p-functions, localized on 
the same site. Nevertheless, their estimate I ~ 1A was 
directly or indirectly used in more later papers . 14 ' 15 ^ 16 In 
fact this integral is estimated to be / « RdpSdpn, where 
Rd P is a cation-anion separation, Sd P -n rfp7r-overlap inte- 
gral. Thus the actual electric polarization induced by the 
spin current is one-two orders of magnitude smaller than 
the authors estimations. 

In Fig[2]we demonstrate the results of numerical calcu- 
lations of several two-site dipole matrix elements against 
3d metal - oxygen separation RmcO- For illustration 
we choose both relatively large integrals (3d z 2\z\2p z ) 
governed by the Me3d-02p cr-bond and the relatively 
small ones (3d xz \z\2p x ) and (3d xz \x\2p z ) governed by 
the Me3d-02p 7r-bond. We make use of hydrogen-like 
radial wave functions with the Clementi-Raimondi ef- 



„ 44.45 ^ e // 



=4.45 and Z 



eff 



=10.53. It is 



fective charges _ ~ 02p — a ^ ^ Me3d - 

clearly seen that given typical cation-anion separations 
RMeO ~ 4 a.u. we arrive at values less than 0.1 a.u. 
even for the largest two-site integral. Reasonable esti- 
mate for the 7r bond integral from the paper by Katsura 
et aL— should be |/(R.dp) | ~ 0.01 A that is two orders of 
magnitude less than that of the authors. 

Relation between local and nonlocal contributions to 
electric polarization is believed to determined by that of 
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^Me3d z2 |z|02p z > 
.<Me3dJx|02p z > 
<Me3d xz |z|02p x > 




FIG. 2: Two-site dipole matrix elements against Me3d-02p 
separation. The arrow near 4 a.u. points to typical Me-O 
separations. 



covalent and overlap effects. The local contribution is 
defined by pure covalent effects and prevails for large 
covalency, that is for large t p d and small E p d, when 
\t P d/E p d\ > S p d- Neglecting the overlap effects we make 
the reliable estimates of nonlocal terms quite question- 
able. 

Interestingly, the nonlocal, or overlap effects are usu- 
ally missed in current calculations of electro-dipole tran- 
sitions in 3d oxides, where one considers the electro- 
magnetic field couples to the electrons via the standard 
Peierls phase transformation of the transfer integral: 



U 



f.. p i(**-*0 



1 he 



A(r)dr, 



(23) 



(24) 



where A is the vector-potential, and integration runs over 
line binding the i and j sites (see, e.g.Reflirf). 

The effective electric polarization differs for the singlet 
and triplet pairing due to a respective singlet-triplet dif- 
ference in the hybridization amplitudes C{ n }(ST). Hence 
we may introduce an effective nonrelativistic exchange- 
dipole spin operator 



p s = n(§x -s 2 ) 

with an exchange- dipole moment 

n = (p> t - cp) s 



(25) 



(26) 



which can be easily deduced from Exps. (TlT)]) and ([2TJ)l . 
For instance, the local contribution of purely oxygen 020 
configuration is 



n 



local 



(Roi + R02) 




(27) 



where E s (p 2 xy ) = 2e x , y + F + ±F 2 , E s (p x p y ) = e x + 
( y + F + ±F 2 , E t (p x p y ) = e x + e v + F - \F 2 are the 
energies of the oxygen two-hole singlet (s) and triplet (t) 
configurations p x ,p 2 and p x Py, respectively, Fq and F 2 
are Slater integrals. We see that this vector is directed 
along y-axis regardless the 5\_ 2 angles and the resultant 
value depends strongly on the CU1-O-CU2 bond geometry 
and crystal field effects. The latter determines the single 
hole energies both for O 2p- and Cu 3d-holes such as 
£xy,xz and e XtV , which values are usually of the order of 
1 eV and 1-3 eV^I respectively. Given estimations for 
different parameters typical for cuprates^ (tpda ~ 1-5 
eV, Fq — 5 eV, F 2 — 6 eV) we can estimate a maximal 
value of ntf 2 c al \ y as 0.01|e|l(~ l0 3 fiC/m 2 ). The local 
contributions to exchange-dipole moment seem to exceed 
the nonlocal ones which are estimated as follows: 



n 



\e\ t -^4(2 Px 



'E 2 (pd) 



x|3eL 



0.001|e|A. (28) 



It is worth noting that for the collinear CU1-O-CU2 bond- 
ing both contributions vanish. As a whole, the exchange- 
dipole moment vanishes, if the M1-O-M2 cluster has a 
center of symmetry. 

Concluding the section it is worth to remind we ad- 
dressed only the charge density redistribution effects for 
Cu 3d and O 2p states, and neglect a direct electronic 
polarization effects for the both metal and anion ions. 
These effects can be incorporated to the theory, if other 
orbitals, e.g. ns- for oxygen ion, will be included to the 
initial orbital basis set. Alternative approach may be 
applied to proceed with these effects, if we turn to a gen- 
eralized shell model42 



IV. RELATIVISTIC MECHANISM OF 
SPIN-DEPENDENT ELECTRIC POLARIZATION 

At variance with a scenario by Katsura et a/3 we have 
applied a conventional procedure to derive an effective 
spin-operator for a relativistic contribution to the elec- 
tric dipole moment in the three-site M1-O-M2 system 
like a technique suggested in refs i 36 ' 37 to derive expres- 
sions for the Cu and O spin-orbital contributions to the 
Dzyaloshinsky-Moriya coupling in cuprates. 

The spin-orbital coupling Vso f° r copper and oxy- 
gen ions drives the singlet-triplet mixing which gives rise 
to a relativistic contribution to electric polarization de- 
duced from an effective spin operator, or an exchange- 
relativistic- dipole moment 



1 



nt = n[§i x s 2 ] , 



(29) 



r sin 2 


(i- 




I 8 







E ?(PxPy) E 2 (p x p y ) 



where Iljj = — i^ool-Pil^ij) IS an exchange-relativistic- 
dipole tensor. It is easy to see that this quantity has a 
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clear physical meaning to be in fact a dipole matrix ele- 
ment for a singlet-triplet electro-dipole transition in our 
three-site cluster . 50 ' 51 ! 52 First of all we should take into 
account the singlet-triplet mixing effects for the ground 
state manifold which are governed by Dzyaloshinsky- 
Moriya interactions 



*s = $s+2j(D-*t); 



T = 



2J 



(30) 



where we make use of Cartesian vector to denote the spin 

triplet function. Then the components of the n tensor 
can be found by projecting P on the spin states 

IT, - -i(9 S \Pi\*Ti) = «*s|i*|*s) - (*T\Pi\* T )) ^J, 

(31) 

In other words, we arrive at a simple form of exchange- 
relativistic-dipole moment as 



p = -in(D- [s x x s 2 ]) 



(32) 



It is worth noting that this vector lies in Cui-0-Cu2 
plane and its direction does not depend on spin config- 
uration. The singlet-triplet overlap density ^>* s *Tj in 
matrix element {* s\Pi\*Tj) has maxima at the points 
R-1,2,3) where the spin-orbital coupling is localized. It 
means that we may pick up the leading local term in 
@2) 



ylocal 



[si x s 2 ]) 



(33) 



((*1*01;10|-P;|*101;10) - (*101;00 \Pi |*1*01;00» , (35) 

where 1*01 labels 101 configuration with d x 2_ y 2 hole on 
Cui site replaced by d* hole with the energy e^* ■ Inter- 
estingly the dipole matrix elements in brackets determine 
the transition probabilities for electro-dipole transition 
d x 2_ y 2 — ► d* on Cui site induced by the covalent and ex- 
change effects in the three-site cluster. Their difference 
can be related with a so called exchange-dipole transition 
moment^ 



P(d->d*) = U(d^d*)(s 1 -s 2 ). 



(36) 



introduced firstly by Y. Tanabe, T. Moriya, and S. 
Sugano 5 - to explain the magnon side bands in 3d mag- 
netic insulators: 



«*1*01;10|P|*101;10) " (*101;Oo|P|*l*01;00» 



n(<4 



(37) 



Interestingly the local contribution to the exchange- 
dipole transition moment vanishes due to the orthogo- 
nality conditions, whereas the nonlocal effects give rise 
both to the in-plane and out-of-plane components both 
of this vector and of the net relativistic electric polariza- 
tion. Indeed, the nonlocal contribution of d x 2_ y 2 — > d* 

spin-orbital excitations on Cui site to the n tensor can 
be written as folows: 



,6 



( d x 2 -v 2 \ l j\df}} ^ 



e 



where II„ and D„ are local (Cui, 2 ,0) contributions to 
the exchange-dipole moment II and Dzyaloshinsky vec- 
tor D, respectively. For a rough estimate we may use 
a relation D/J ~ Ag/g, where g is the gyromagnetic 
ratio and Ag is its deviation from the value for a free 
electron^ 

Another contribution to IT, = — i{*s\Pi\*Tj) we ob- 
tain, if make use of singlet and triplet hybrid functions 
^ioi;SM perturbed by spin-orbital coupling as follows: 36 



^Wl-.SM = ^10X;SM — 



E 



(*M;r<S'Af |^o|*101: 



SM 



\*}&M'T' E*s> + i T ,({n}) - S as+ ir o (101) 



{n};T'S'M' 



(34) 

Notice that {n} for the hybrid function ^{„} ; r'S'M' 
points only to a bare, or generative ionic configuration. 

For an illustration we address the 3d x 2_ y 2 — > 3<i* exci- 
tations driven by Vso{Cu\) within ground state 101 con- 
figuration. The proper contribution to the singlet-triplet 
matrix element of P can be written as follows 



n, 



101:00 



PA* 



101 



{d a 



\lj\d*) 



E t {dp a ) - ep E s {dp a ) - ep 



{2p a \xi\3d 



(i) 



> + 



E t {dp a ) E s 



^-yj (2p a \x i \3d§ ) )] , (38) 



thus we arrive at nonzero H xz , Hyz components provided 
d* = d xy and U zy component provided d* — d xz , if to 
account for the nonvanishing overlap dipole matrix ele- 
ments (2p a \x a \3d x 2_ y 2) and (2p x \z\3d xz ). A reasonable 
estimate for the maximal value of IT, can be made, if 
address the relation (J55J): |n y | - 0.1II - 10~ 4 |e|A. 

It should be noted that for the contribution of bare 
configurations other than that of ground state 101 we 
may use a simplified expression 36 



^101;SM = &101:SM + 



{n}r 



C{n}{ 



2S+1 T 



{n}\TSM 



E 



(4> 



{n};T'S'M' \Vso\^ {n} -r S m) 



s fa r , E 2S , +lr ,{{n}) - E 2S+1To {101) 



{n}-T'S'M' 



(39) 
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However, on closer examination we arrive at vanish- 
ing contribution of these terms to exchange-relativistic- 
dipole moment. 

Thus the Dzyaloshinsky-Moriya type exchange- 
relativistic-dipole moment (|32|) is believed to be a dom- 
inant relativistic contribution to electric polarization in 
Cui-0-Cu2 cluster. It is worth noting that the exchange- 
dipole moment operator (|25p and exchange-relativistic- 
dipole moment operator (|32|) are obvious counterparts of 
the Heisenberg symmetric exchange and Dzyaloshinsky- 
Moriya antisymmetric exchange, respectively. Hence, the 
Moriya like relation \Hij\ ~ Ag/g\U\ seems to be a rea- 
sonable estimation for the resultant relativistic contri- 
bution to electric polarization in M1-O-M2 clusters. At 
present, it is a difficult and, probably, hopeless task to 
propose a more reliable and so physically clear estima- 
tion. Taking into account the typical value of Ag/g ^0.1 
we can estimate the maximal value of |IIy | as 10 -3 |e|A(~ 
10 2 /iC/m 2 ) that points to the exchange-relativistic mech- 
anism to be a weak contributor to a giant multiferroicity 
with ferroelectric polarization of the order of 10 3 /xC/m 2 
as in TbMnOa, 1 though it may be a noticeable contrib- 
utor in e.g. Ni 3 V 2 8 3 

V. PARITY BREAKING EXCHANGE 
COUPLING AND EXCHANGE-INDUCED 
ELECTRIC POLARIZATION 

Along with many advantages of the three-site cluster 
model it has a clear imperfection not uncovering a direct 
role played by exchange coupling as a driving force to in- 
duce a spin-dependent electric polarization. Below we'll 
address an alternative approach starting with a spin cen- 
ter such as MeO„ cluster in 3d oxides exchange-coupled 
with a magnetic surroundings. Then the magnetoelectric 
coupling can be related to the spin-dependent electric 
fields generated by a spin surroundings in a magnetic 
crystal. In this connection we should point out some 
properties of exchange interaction that usually are missed 
in conventional treatment of Heisenberg exchange cou- 
pling^ Following after paper by Tanabe el al£ (see, also 
Refll2l) we start with a simple introduction to exchange- 
induced electric polarization effects. 

Let address the one-particle (electron/hole) center in a 
crystallographically centrosymmetric position of a mag- 
netic crystal. Then all the particle states can be of defi- 
nite spatial parity, even (g) or odd (u), respectively. Hav- 
ing in mind the 3d centers we'll assume the even-parity 
ground state \g). For simplicity we restrict ourselves by 
only one excited odd-parity state \u). The exchange cou- 
pling with surrounding spins can be written as follows: 

VU = 5^J(R n )(s-S B ) ) (40) 

n 

where J(R n ) is an orbital operator with a matrix 

i<^) = ( I T Bg( £r\ HJHY ( 41 ) 



The crystallographic centrosymmetry condition requires 
that 

^V(R„) = ]T/ n3 (R n ) = 0. (42) 

n n 

The parity-breaking off-diagonal part of exchange cou- 
pling can lift the center of symmetry and mix \g) and \u) 
states 

\g) ^ \g) + c gu \u) , (43) 

where 

Csu = A- g 1 ^/ ffU (R„)(s-S„) (44) 

n 

with A ug = e u — e g . In turn, it results in a nonzero 
electric dipole polarization of the ground state 

P = 2c 9 „( 9 |er|t t )=^n„(s-S n ), (45) 

n 

where d = er is a dipole moment operator, 

n„ = 2I gv [K n )—r . (46) 

It is easy to see that in frames of a mean-field approxima- 
tion the nonzero dipole moment shows up only for spin- 
noncentrosymmetric surrounding, that is if the condition 
(S(R„)) = (S(— R„)) is broken. For isotropic bilinear 
exchange coupling this implies a spin frustration. 

Kinetic contributions to conventional diagonal and un- 
conventional off-diagonal exchange integrals can be ob- 
tained, if one assume that surrounding spins are formed 
by a single electron localized in the same \g) state: 




where t gg is a transfer integral between ground \g) states 
of the neighboring ions, while t ug is a transfer integral 
between ground \g) state of the neighboring ion and \u) 
state of the central ion, A g9 is the energy of the charge 
transfer between ground \g) states of the neighboring 
ions. 

It should be noted that at variance with DM type 
mechanism the direction of the exchange-induced dipole 
moment for i,j pair does not depend on the direction 
of spins Si and Sj. In other words, the spin-correlation 
factor (Si • Sj) modulates a pre-existing dipole moment 
II which direction and value depend on the Mej-O-Mej 
bond geometry and orbitals involved in exchange cou- 
pling. 
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The net exchange induced polarization of the magnetic 
crystal depends both on crystal symmetry and spin struc- 
ture. The allowed direction of the average P in crystal 
can be unambiguously determined by symmetry analysis, 
for instance, P should be parallel to all the mirror planes 
and glide planes. 

The magnitude of off-diagonal exchange integrals can 
sufficiently exceed that of conventional diagonal exchange 
integral mainly due to a smaller value of the energy sepa- 
ration A gg — A ug as compared with A gg and larger value 
of transfer integral t ug as compared with t gg due to the 
purely oxygen character of odd-parity \u) state. Given 
reasonable estimations for off-diagonal exchange integrals 
I ug 0.1 eV, g — u energy separation A gu sa 2 eV, 
dipole matrix element |(<?|er|u)| ss O.lA, spin function 
|((s • S„))| 1 we arrive at estimation of maximal value 
of electric polarization: P w 10 4 /zC/m 2 . This estimate 
points to exchange-induced electric polarization to be a 
potentially the most significant source of magnetoelectric 
coupling for new giant multiferroics. 

It is worth noting that the exchange-induced polar- 
ization effect we consider is particularly strong for the 
3d clusters such as MeO„ with the intensive low-lying 
electro-dipole allowed transition \g) — ► \u) which both 
initial and final states are coupled due to a strong ex- 
change interaction with a spin surroundings. This sim- 
ple rule may be practically used to seek new multiferroic 
materials. 

The parity-breaking exchange coupling can produce a 
strong electric polarization of oxygen ions in 3d oxides 
which can be written as follows 

Po=^II n «S >-S n ), (49) 

n 

where S„ are spins of surrounding 3d ions, (So) oc 
In^ra is a spin polarization of oxygen ion due to sur- 
rounding 3d ions with I n being the exchange coupling 
tensor. It seems the oxygen exchange-induced electric 
polarization of purely electron origin is too little appreci- 
ated in the current pictures of multiferroicity in 3d oxides. 

VI. LACK OF SPIN-DEPENDENT ELECTRIC 
POLARIZATION IN EDGE-SHARING Cu0 2 
CHAINS 

According to the phenomenological theory by 
Mostovoyi - and microscopic model by Katsura et al£ 
the spin-spiral chain cuprates LiCuV04 and LiCu2 02 
seem to be prototypical examples of ID spiral-magnetic 
ferroelectrics revealing the relativistic mechanism of 
" ferroelectricity caused by spin-currents". Indeed, 
the net nonrelativistic polarization of a spin chain 
formed by Me 3d ions even with no center of symmetry 
inbetween can be written as follows 5 

p e// = n Y, • s i+ i) - (s, • s^)] , (so) 

j=even 




FIG. 3: The fragment of a typical edge-shared Cu02 chain. 
Note the antiparallel orientation of the oxygen Dzyaloshinsky 
vectors directed perpendicular to the chain plane. 



hence for a simple plane spiral ordering both the on- 
site and net polarizations vanish while the spin-current 
mechanism-'' 10 directly points to a nonzero polarization 
concomitant spin spiral order. However, a detailed anal- 
ysis of relativistic effects for the system of e g -holes in 
a perfect chain structure of edge-shared Cu04 plaque- 
ttes as in LiCuV04 shows that the in-chain spin cur- 
rent does not produce an electric polarization. First of 
all we should point to a high symmetry of CU1-O-CU2 
bonds in edge-sharing Cu02 chains (see Fig. [3]) that re- 
sults in a full cancellation of a net Dzyaloshinsky vector, 
though the partial oxygen contributions survive being of 
opposite sense . 36 i 37 Cancellation of the Dzyaloshinsky- 
Moriya coupling in perfect edge-sharing Cu02 chains im- 
plies immediately the same effect for the net exchange- 
relativistic-dipole moment P. Indeed, the dominant con- 
tribution to the exchange-relativistic-dipole moment for 
isolated CU1-O-CU2 bonds is governed straightforwardly 
by the respective Dzyaloshinsky vectors, hence their can- 
cellation for CU1-O/-CU2 and Cui-0//-Cu2 bonds in 
edge-sharing Cu02 chain geometry (see Fig|3]) leads to 
the vanishing of the exchange-relativistic electric polar- 
ization. It seems, small nonlocal terms addressed in Sec. 
IV could survive, however, the symmetry considerations 
point to their vanishing as well. Indeed, both the xz and 
yz components of the Ily tensor differ in sign for the 
CU1-O/-CU2 and CU1-O//-CU2 bonds while the zy com- 
ponents differ in sign for the contribution of Vso(Cu\) 
and Vso{Cu2)- Thus we may state that the edge-shared 
CUO4 plaquettes chain arrangement appears to be ro- 
bust regarding the proper spin-induced electric polariza- 
tion both of the nonrelativistic and relativistic origin. It 
means that we should look for the origin of puzzling mul- 
tiferroicity observed in LiCuV04 and LiCu202 some- 
where within the out-of-chain stuff. 54,55 
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VII. CONCLUSION 

We have considered a microscopic theory of spin- 
dependent electric polarization in 3d oxides starting with 
a generic three-site two-hole cluster. A perturbation 
scheme realistic for 3d oxides is applied which implies the 
quenching of orbital moments by low-symmetry crystal 
field, strong intra-atomic correlations, the pd-transfer ef- 
fects, and rather small spin-orbital coupling. An effective 
spin operator of the electric dipole moment is deduced in- 
corporating both nonrelativistic oc (§i-§2) and relativistic 
cx [§i x §2] terms. The nonrelativistic exchange-dipole 
moment is mainly governed by the effects of the redistri- 
bution of the local on-site charge density due to pd cova- 
lency and exchange coupling. The relativistic exchange- 
dipole moment is mainly stems from the nonrelativis- 
tic one due to the perturbation effect of Dzyaloshinsky- 
Moriya coupling and is estimated to be a weak contrib- 
utor to the electric polarization observed in the most 
of 3d multiferroics. Our description is focused on Cui- 
O-CU2 clusters typical for different cuprates, however, 
the generalization of the results onto the M1-O-M2 clus- 
ters in other 3d oxides is trivial. The approach real- 
ized in the paper has much in common with the mech- 
anism of the bond- and site-centered charge order com- 
petition (see, e.g. RefHa) though we start with the el- 
ementary pd charge transfer rather than the dd charge 
transfer. An alternative approach to the derivation of the 



spin-dependent electric polarization is considered which 
is based on the parity-breaking exchange coupling and 
exchange induced polarization. 

We point to the oxygen electric polarization effects due 
to an exchange-induced electric fields to be an important 
participant of the multiferroic performance. Anycase, the 
nonrelativistic electronic polarization mechanism is be- 
lieved to govern the multiferroic behaviour in 3d oxides. 

It is shown that the perfect chain structure of edge- 
shared CuC>4 plaquettes as in LiCuV04 or LiCu2C>2 ap- 
pears to be robust regarding the proper spin-induced 
electric polarization both of nonrelativistic and relativis- 
tic origin. In other words, in contrast with the predic- 
tions of the model by Katsura et al£ the in-chain spin 
current does not produce an electric polarization. Hence 
the puzzling multiferroicity observed in LiCuV04 and 
LiCri202 18,19 originates from an out-of-chain stuff. 

Clearly, the model approach applied can provide only 
a semiquantitative description of magnetoelectric effects 
in 3d oxides. More correct account for the overlap, or 
nonorthogonality effects and those produced by a non- 
magnetic surroundings of the three-site two-hole cluster 
are needed. 
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